Criteria for dynamically stable decoherence-free subspaces and incoherently generated 

coherences 
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We present a detailed analysis of decoherence free subspaces and develop a rigorous theory that 
provides necessary and sufficient conditions for dynamically stable decoherence free subspaces. This 
allows us to identify a special class of decoherence free states which rely on incoherent generation 
of coherences. We provide examples of physical systems that support such states. Our approach 
employs Markovian master equations and applies primarily to finite-dimensional quantum systems. 
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I. INTRODUCTION 

The theory of decoherence-free (or noiseless) subspaces 
(DFS) provides an important strategy for the passive 
preservation of quantum information [J, 0, and has 
been subjected to successful but restricted experimental 
tests H S S B i, DFS theory has been developed 
in the context of algebraic symmetries in the interaction 
Hamiltonian by Zanardi and Rasetti (ZR) 0, [ll| and by 
Lidar, Bacon, and Whaley and in the context of 
semigroup dynamics via quantum master equations by 
Lidar, Ch uang , and Whaley (LCW) Q and by Shabani 
and Lidar [ij]. The latter analysis has been extended to 
the theory of DFS in the non-Markovian regime [ij, [isj 
and to DFS with imperfect initialization |l4|. Analogues 
of DFS have also been proposed in the quantum channel 
formulation [l^. Potential advantages of DFS include 
quantum circuit simplifications by reducing the demands 
for quantum error correction and for quantum memory 
and transport. Given the significance of DFS for realiza- 
tion of quantum information processing, it is important 
that the foundations of the theory are rigorous and unam- 
biguous, especially in preparation of future experiments 
to test and exploit the DFS. 

DFS have been defined as collections of states that 
undergo unitary evolution in the presence of couplings 
to the environment (i.e., decoherence effects). However, 
unitary evolution of a quantum state can arise in a num- 
ber of ways and this fact has resulted in the development 
of related, but different definitions of DFS in the litera- 
ture. In the context of Markovian master equations, DFS 
have frequently been defined as a collection of states for 
which dissipative (decoherence) part of the Markovian 
master equation is zero. In this paper, we give a detailed 
analysis of DFS based on Markovian semigroup master 
equations defined for finite-dimensional Hilbert spaces. 
We provide a definition for dynamically stable DFS that 
identifies a new subclass of DFS states. This definition 
allows us to develop a theorem that gives necessary and 
sufficient conditions for the existence of dynamically sta- 



ble DFS and that provides a constructive protocol for 
determining the composition of all DFS. 

Existence of this new class of DFS states is based on 
the rather counter-intuitive phenomenon that the deco- 
herence term, which is of Lindblad type, can generate 
coherences between these states and other states which 
are then exactly cancelled by the system Hamiltonian. 
This is in striking contradiction to the usual role of the 
decoherence term, which includes damping of coherences 
rather than cancellation of these. This new phenomenon 
is different from techniques to suppress decoherence like 
dynamical decou pling (itI . [isl . [l9l . [20| and continuous er- 
ror correction [2l|, [22l because it does not rely on the rel- 
ative time scale for the system evolution and decoherence 
processes. We show that the existence of these incoher- 
ently generated coherences results from an invariance of 
the Markovian master equation under a specific trans- 
formation and provide two examples of physical systems 
that support these special decoherence-free (DF) states. 
We also indicate that this new phenomenon is rare and 
does not exist for finite-dimensional systems subject to 
the common quantum noise processes such as sponta- 
neous emission, thermal excitation and dephasing. 



II. BACKGROUND 

Our analysis of the DFS is predicated on treating the 
open system as Markovian and thus using Markovian 
master equations. 

A system S with state vectors in finite-dimensional 
Hilbert space Tis is coupled to a reservoir R with state 
vectors in Hilbert space Hr. The system-plus-reservoir 
(S-l-R) dynamics are fully described by the Hamiltonian 
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(2.1) 



with Ha,, Hy{ and iJjnt the system, reservoir and inter- 
action Hamiltonians, respectively, and 1 the identity op- 
erator. We are interested in the evolution of the sys- 
tem's reduced state p, which is obtained by tracing over 
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the reservoir's degrees of freedom (subject to appropri- 
ate assumptions: no initial correlations between S and a 
Markovian R and weak SR coupling). Dynamics of the 
reduced density matrix p depends on the effective-system 
Hamiltonian H^s = + A, with A the Hermitian con- 
tribution that stems from the interaction between the 
quantum system and the reservoir. The evolution of the 
reduced density matrix p also depends on the dccoher- 
ence super-operator 



k.l=l 



aki 



(2.2) 



where A = {au) is a positive semidefinite and time- 
independent matrix and No, — dim(7is)- The set of op- 
erators 



} , Fo = i, (2.3) 



is a basis for the space of all linear transformations de- 
fined on Tis- 

The evolution of p is given by the master equation 



P = -i[Hcii,p] + Ld[p]^ 



(2.4) 



Henceforth, we refer to Eq. (|2.4|) as Markovian open sys - 
tem dynamics. Eq. (|2.4p is of Lindblad form |23l . |24| . 
This means that Eq. (|2.4p generates a completely-positive 
trace preserving map and thereby ensures that the solu- 
tion to the master equation has all the required properties 
of a physical density matrix at all times. 

We note here that the terms A and Ld [p] in Eq. (|2.4[) 
originate from the system-reservoir coupling and induce 
a modified unitary evolution and decoherence in system 
S. In this paper we are interested in characterizing states 
that undergo purely unitary time evolution. 



III. RESTRICTED DECOHERENCE-FREE 
SUBSPACES 

DFS have been defined as collections of states that 
evolve in a unitary fashion 0, i, i, [Til, [H, [H, [13. We 
shall refer to the complete set of such states as a dynam- 
ically stable DFS. From now on we will always assume 
that the state of the system p{t) satisfies Eq. (|2.4p . In 
this section, we characterize the restricted subset of DFS 
states p fulfilling Lj:,[p{t)] = for all t. Conditions for the 
more general definition of all states evolving in a unitary 
fashion will be considered in Sec. IIVI 

The traditional approach to DFS in the context of 
Markovian master equations requires that Lj:,[p(t)] = 
This constitutes a condition for an instantaneous 
DFS, i.e., a DFS at a specific time t. However, satisfying 
the condition LY)[p{t)] = at one specific time alone, i.e., 
having an instantaneous DFS, does not guarantee unitary 
evolution. The action of the effective system Hamilto- 
nian Hcff modifies the density matrix p{t) and thus at a 



later time t' > t, p{t') might no longer describe a pure 
decoherence-free state. The condition L-£,[p{t)] = for all 
t further ensures that state p{t) undergoes pure unitary 
evolution and thus extends the notion of an instantaneous 
DFS. Nevertheless, as we will demonstrate explicitly in 
Section 13 this condition is sufficient, but not necessary, 
for the existence of a DFS. Therefore we shall refer to 
the set of states satisfying this sufficiency condition as a 
restricted DFS. 

Definition 1. Let the time evolution of a state p{t) be 
given by Markovian open system dynamics. Then a state 
p{t) is a pure restricted decoherence-free state if p^it) = 
p{t) and LdW)] = 0, \/t. 

Remark 1. We require states to be pure in Definition [1] 
i.e., p = IV') ("01 a-nd p^{t) = p{t), in order to ensure 
that all states satisfying this definition do indeed un- 
dergo unitary evolution. For mixed states, it is possi- 
ble to find states that are formally decoherence-free (i.e., 
Lnipit)] = for all t), yet do not undergo a proper uni- 
tary dynamics. Examples of such states include the com- 
pletely mixed state p{t) = Is/A^s in a system that deco- 
heres through dephasing only, as well as many thermal 
states. The reason that these states are not decohering 
according to Definition [1] is that the incoherent transition 
rates between any two states in the mixture are equal, 
without the dynamics of each being unitary. Such states 
contain minimal information and are not useful for the 
purpose of quantum information processing. Hence we 
focus in this paper on pure states. Our results can be 
generalized to specific mixtures of pure decoherence-free 
states 3 but not to a thermal state. 

This definition of a restricted DFS implies that a pure 
state p{t) would remain pure during its evolution and 
thus be an ideal tool to process quantum information. It 
is therefore important to find criteria for the existence of 
a restricted DFS in a given quantum system. We pro- 
vide such criteria in Section [ill Bl below. Before deriving 
these criteria, we first discuss the conditions for an in- 
stantaneous DFS for which LD[p(i)] = at a fixed time 
t, tightening existing conditions that have been derived 
previously in the Markovian setting Q . 



A. Instantaneous decoherence-free states 

The condition for existence of instantaneous DFS 
states has been previously studied by LCW, who formu- 
lated the following theorem which we restate here and 
then tighten. 

Theorem 1. (LCW Theorem 0]) A necessary and 
sufficient condition for generic decoherence-free dynam- 
ics (i.e. LoIp] =: 0) in a subspace {TL — span[{\i)W of 
the register Hilbert space is that all basis states \i) are 
degenerate eigenstates of all the error generators T : 



Fk\i) = Ck\i) V fc, i. 



(3.1) 



3 



// the error generators {Fk} can be closed as a semi- 
simple Lie algebra, then all states in Ti. are annihilated 
by all error generators T , i.e. Fy;\i) = 0|z) V /c, i. 

We first identify and remedy shortcomings correspond- 
ing to the fohowing two examples concerning: i) zero 
eigenvalues for the matrix [au) in Eq. (|2.2p : and ii) non- 
zero eigenvalues for the error generators in Eq. (|3.ip . 
These two examples demonstrate the associated difficul- 
ties with this description. 

Example 1. Consider a 3-level system with basis 
{|0), |1), 1 2)} with vanishing Hamiltonian Hgff and Ld[p\ 
of the form 112. 2\) with Ng = 3. A is an 8 x 8 matrix 
because there are iV| — 1 = 8 error generators in this 
example. We choose the components of A as Oki = 1 
for k,l = 1,2 and zero otherwise, so that only the error 
generators Fi = |0)(1| and F2 = |0)(2| are relevant. We 
can rewrite A as 



with 



A = BKB\ B 
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A= I ^ 
O Oj 
(3.2) 

where O denotes sub-matrices with all components equal 
to zero and 1 is a 6 x 6 identity matrix. Obviously A 
possesses zero eigenvalues. The state = 2^"'^/^(|l) — 
1 2)) is in a DFS because 



LnU){i^\] = [H,ff,\^){i^\]=Q. 



(3.3) 



On the other hand we have Fi^2\^) = ±2 "'^/^jO) 7^ 
ci,2\tp) , SO that \ip) is not an eigenstate of the error gen- 
erators Fk, and the LOW theorem would thus not identify 
\tp} as a DFS state. 

This problem originates from the inclusion of zero 
eigenvalues of A [401 '^^^ be resolved as follows. 
Let B = (bki) be a unitary matrix that diagonalizes the 
matrix A and let A^- > be the positive eigenvalues of 
matrix A. Then we can define Lindblad operators 



(3.4) 



fe=i 



Ji 



F1+F2 
V2 



|0)(1| + |0)(2| 
V2 ■ 



(3.7) 



Then Ji\'ip) = 0|V'). Thus the LCW theorem holds if 
it is modified to require that \^) is an eigenstate of the 
operators {J/} rather than of {Fk}. 

Example 2. Consider a 2-level system with basis 
{|0),|1)} and Ld[p] of the form with M = 2, 

Ai = A2 = 1 and the decoherence operators 



ii = |i)(i|,i2 = |o)(o| + |o)(i|. 



(3.8) 



In order to close the Lie algebra C generated by operators 
Ji and J2, we need to introduce the operator J3 = |0)(1|. 
The sub-algebra spanned by J3 forms an Abelian ideal so 
that C is not a semi-simple Lie algebra. The state |0) is 
an eigenstate of all decoherence operators, Ji\0) = 0\0), 
J2IO) = 1|0), J3IO) = 0|0). Therefore, according to the 
LCW Theorem, |0) should be decoherence- free. However, 



Ld[\0){0\] 



|0)(1| + |1)(0| 



^0. 



(3.9) 



Hence |0) does not exhibit decoherence- free dynamics, yet 
it is an eigenstate of all the error generators. 

This example shows that, in addition to excluding zero 
eigenvalues of A, another condition for eigenstates of er- 
ror generators Ji corresponding to non-zero eigenvalues 
is required in order to tighten the LCW theorem. Here 
we establish a simple and complete criterion for identi- 
fying whether a given state is an instantaneous DFS for 
a general decoherence process described by a Markovian 
master equation. We begin with a definition of the de- 
coherence operator F for a system having decoherence 
super-operator ^dH, Eq. (|3.5p . of the form 



M 



t^Y.^iJlJi. 



(3.10) 



1=1 



The eigenvalues of f correspond to the inverse lifetimes 
of the corresponding eigenstates. 



that transform Eq. 
ence super-operator 



to the 'diagonalized' decoher- Proposition 2. 



Lt>[p] 



M 

-Y 

2 ^ 

(=1 



LDU)m = ^ 



(3.11) 



h ({Jupj}] + [JipJ}]) , (3.5) iff M-^) = qIV') for all I = 1,...,M and V\^) = 5]^), 

^ ^ , > v^M \ I ,2 



where we have relabeled the eigenvalues {A;} and opera- 
tors {J;} so that \i > for / = 1, . . . , M < A^l - 1 and 
A; = otherwise. In this form the zero eigenvalues of A 
do not contribute to ioM- 

In the above example we then have 



where g = Y^i^i Mci 

Proof. First we prove that the conditions in Proposition 
[5] are sufficient. Suppose there exists a state such 
that Ji\i}) = ci\ip) and f|V') = .g|V')- Then 



L^[p]^2{[J^, pJl] + [J,p,Jl]) 



(3.6) 



M / M \ 

^A,j,i^)(V'iJ/ = ^A,iQn 
1=1 \i=i ) 



Mm. (3.12) 
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and for p = \ip) {ipl 



(3.13) 



because g = Y^fii ^^/I^P- 

To prove that the conditions in the proposition are 
necessary, we now assume Lb[p] = 0. First we show that 
Jilip) = ci\ip). Using Eq. (|3.5p wc evaluate the expression 



M 



1=1 



(3.14) 

Jilijj) can generally be written as Ji\ip) = ci\ip) + 
with \ipi~) some (non- normalized) state that is orthogonal 
to the state Substituting this into Eq. (|3.14p yields 



M 



(3.15) 



1=1 



Because A; > for ^ = 1, . . . , M we find || |V'J-)i|l ~ Oj i-^^-j 
= ci\ip) I = 1, . . . , M . We also need to show that 
= .9|V')- Consider a subspace 

|(V'±IV^)=0}. (3.16) 

Then 

= {4>i_\LT,[pm = -hi,^\t\i^) V|Vi) e V^. (3.17) 



This is only possible if f IV") = fflV')- 

The value of g can be derived as follows. Consider the 
operator 



M 



(3.18) 



1=1 



The state jV') is an eigenstatc of operators G and with 
eigenvalues g' = X]f=i -^ik'P ^^.'^ 5i respectively. Then 



(5')' = I|G^K')IP = (^|GG^I^>=53', 



(3.19) 



assuming that the state \^) is normalized. Therefore, 



9 = TaLiMci 



□ 



Remark 2. Our Proposition [2] immediately removes ex- 
ample [2] because the example does not satisfy the condi- 
tion flV-) = sl^)- 

Corollary 3. Suppose that all eigenvalues of operators 
Ji for I = 1, . . . ,M are equal to zero; then the state jf/') is 
an instantaneous pure decoherence-free state iff it is an 
eigenvector of the decoherence operator \3.10\) with an 
eigenvalue of zero. 



We now discuss two examples for which the second con- 
dition in Proposition [21 rji/)) = g\'ip), is rendered redun- 
dant by virtue of the intrinsic physics, i.e., the physical 
system automatically satisfies this condition. These ex- 
amples are drawn from the most common types of noise 
in quantum systems, namely excitation, decay and de- 
phasing. 



The first example is noise due to excitation or decay. 
During excitation, the system transits to a different state 
and gains energy. For decay, the system changes its state 
and looses energy. For decay processes the correspond- 
ing operators Ji can be represented by upper triangu- 
lar matrix with zeros along the diagonal and thus can 
have only zeros for its eigenvalues. Excitation processes 
are correspondingly represented by lower triangular ma- 
trices with zeros along the diagonal and therefore also 
have only zero eigenvalues. Consequently this kind of 
noise automatically satisfies Corollary [3] and hence the 
second condition in Proposition [2] ( i.e., r|'0) = 5|'0)) is 
redundant for quantum systems with noise comprised of 
independent decay and excitation processes. 



The second example where the second condition in 
Proposition [5] becomes redundant is pure dephasing. 
During dephasing [soj , a quantum system remains in the 
same state, but acquires a phase. Thus for quantum noise 
due to dephasing, there exists a representation such that 
all Ji are diagonal. In this case, J/ commutes with J^^, 
which means that J/ is a normal operator [26j . As a re- 
sult, we know that eigcnstates of J; are also eigenstates of 
jj . Consequently the action of each term contribut- 
ing to f in Eq. (PTU)) is of the form j/ Ji|-0) ^qI^/') 
with a simultaneous eigenstate of all J; and j/, so 
that we automatically have T\'ip) = g\'ip) with g a con- 
stant. Thus in this situation the second condition of 
Proposition [2] also becomes redundant. 



At this point it is useful to compare the results of Za- 
nardi [J to the LCW theorem and this new tightened 
version. In Zanardi developed a theorem similar to 
theorem 3 of LCW @. Zanardi's theorem applies only 
to a restricted set of systems and does not need tight- 
ening. The theorem refers to a quantum system S that 
consists of identical spin-like particles, each of which pos- 
sesses a finite number of equally spaced energy levels. It 
is also assumed that this system admits decoherence due 
to either dephasing, or decay, or excitation, alone and 
not to a combination of these processes. As we indicated 
above, the new condition in Proposition [5] is not required 
for such independent sources of decoherence. In general. 
Theorem [1] with Proposition [2] provides tight conditions 
for the existence of instantaneous DFS in all situations 
described by a Markovian master equation. 
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B. 



Existence criteria for dynamically stable 
restricted DFS 



In section UlI A[ we determined when pure states p{t) 
are decoherence-frec at a fixed time t. However, as men- 
tioned before, this is not sufficient for unitary dynamics 
because the effective system Hamiltonian iJoff can drive 
p{t) out of the decoherence-free subspace at a later time 
t. In this section we will incorporate time evolution in 
our considerations in order to establish more precise cri- 
teria for the existence of DFS at all times. This will 
characterize restricted DFS that satisfy L£i[p{t)] = OVt, 
i.e., dynamically stable DFS. The following theorem es- 
tablishes the composition of TYdfs in terms of all initial 
states that undergo decoherence-free evolution. 

Theorem 4. Let the time evolution be given by Marko- 
vian open system dynamics, Eq. \2.4-^ . The space 



span \lp2), . . . , IV'if)} = Ejurap 



(3.20) 



is a DFS for all time t iff Ejump is a subspace that is 
invariant under H^g and that has basis vectors satisfying 



(3.21) 



for all I ~ 1 , . . . , Af and for all k ~ I, . . . , K with g ~ 

Proof. (Note: we have used the notation -Bjump to em- 
phasize the relation with the operators J which corre- 
spond to the jump operators in a quantum trajectory 
description (sec below).) To prove that conditions of 
Theorem 2] are necessary, we suppose i?jump is a DFS. 
Then any state 



IV'(O)) e Hbfs 



can be written as 



K 



|^(0))=5]afe|^fe). 



fc=i 



which evolves in time according to 

\m) = u{t)\m) 



(3.22) 



(3.23) 



(3.24) 



with U{t) — exp(— liifeff) because Ld — 
OVi. As 

mt)} e HdFS = ^Jump, (3.25) 

we can write the evolving state as 

K 

|V'(t))=5]afe(t)|Vfe). (3.26) 

k=l 

As the state \ip{0)) e Hufs is arbitrary, this implies that 



K 



K 



(3.27) 



i.e., the Hamiltonian HcS leaves TiuFS invariant. 

Since LDmO))mO)\] = for aU |^(0)) G Ej^^^, by 
Proposition [2] we know that basis states \tpk) must be 
eigenstates of J/ and T V I = 1, . . . , M. To complete the 
proof that the conditions of Theorem 2] are necessary, 
we have to show that the eigenvalues c; and g are equal 
for all basis states \ipk)- We do this by establishing a 
contradiction. Suppose two arbitrary basis states \ipk) 
and {ipk') have different eigenvalues: 



Jl\ipk) = Q,fe|V'fe)) Jlli^k') = Q,fe'|V'fc')- 
Then the state 

v2 



(3.28) 



(3.29) 



is not an eigenstate of Ji. However, because \ip) G TYdfs 
Proposition [5] implies that must be an eigenstate. 
Hence the eigenvalues must be equal for all basis states. 

To prove that the conditions of Theorem |4] are suffi- 
cient, we consider the action of the Liouvillcan 



C[.] = -i[H,s, .]+Lo[: 



(3.30) 



on a coherence |V'o)(V'ol between two arbitrary states 
\i/jo), \ipo) € £^jump- We then have 



0, 



(3.31) 



because the states are eigenstates of Ji and f . Hence the 
action of C can be reduced to 



-i[H,s, |^o)(V^ol] 

-iiv^i)(v^i+zivo)(^;i 



(3.32) 



with -ffcftlV'o) = iV'i)- Because Hcs leaves i?jump invari- 
ant we know that G i?jump. This is a superpo- 
sition of two coherences that are of the same type as the 
original coherence. Repeated application of £ and the 
linearity of Ld allows us to infer that the action of Ld 
will not contribute to any power £" [p] . Consequently, we 
find for the time evolution of an arbitrary density matrix 
p(0) formed of states G i?,jump, 



[p(0)] = e 



P(0) 



«p(0)e 



(3.33) 



fc=i 



k=l 



Hence, p{t) e Ti-uFS, since it shows unitary evolution 
under HcS- □ 

Remark 3. By adopting Definition [T] for DFS, we were 
able to derive that i?jump is invariant under action of H^g 
as a consequence of the definition, obviating the need to 
impose this condition as part of the definition itself, as 
was done in Ref. [lli. 
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Next we provide a Proposition with conditions for 
^'jump that are equivalent to requiring -Ejump to be invari- 
ant under i?eff, i-e., to one of the conditions in Theorem^] 
above. The advantage of these alternative conditions is 
that it might be easier to verify them for systems that 
support DFS with large dimensions. 

Proposition 5. Let Ej^mp be a space spanned by some 
or all of a collection of simultaneous degenerate eigen- 
states of Ji with eigenvalues ci for I = 1, . . . , M that are 
also eigenstates of T with eigenvalue g. Then Ej^mp is 
invariant under H^ff iff 



Heff, Ji 1^) = IV-) for all 1^1,..., M (3.34) 



pure state, i.e., with Tr[p^(0)] = 1). This motivates us to 
provide a more general definition for the DFS as a collec- 
tion of states that evolves in a unitary fashion. This is a 
larger set than the restricted DFS studied in Section Hill 
Let D{'H) be the set of all density matrices that de- 
scribe pure states and that arc defined for a quantum 
system associated with the Hilbert space H. 

Definition 2. Let the time evolution of an open quan- 
tum system with Hilbert space Tis be given by Markovian 
open system dynamics, Eq. \2.4^ . Then a decoherence 
free subspace TioFS is a subspace ofTis such that all pure 
states p(t) e D{Hdfs) fulfill 

dtTr[p^{t)] = OVt > 0, with Tr[p^{0)] ^ I. (4.1) 



(3.35) 



for all IVj) G Ejump- 



Proof. Let be an eigenstate of j; (/ = 1, . . . , M) and 
r, with respective eigenvalues c; and g. Suppose Hcs 
leaves -Bjump invariant. Then Hes\4') is an eigenstate of 
operators Ji and T with respective eigenvalues ci and g. 
Consequently, Eqs. (|3.34p and (|3.35p hold. Suppose now 
that Eqs. (|3.34|) and (|3.35|) are true for any & i?jump- 
Then 



(3.36) 



t (HcsH)) = g (hM) , 



i.e. Hcs\ip) is in £;jump for all G Ej 



ump • 



□ 



Remark 4. For the specific examples of decoherence pro- 
cesses describing either dephasing, or decay, or excitation 
in a finite-dimensional spin-like systems, it was already 
shown in Ref. Q that the set of states 

{|V'i>,|^2),...,|^k)} (3.37) 

that fulfills 

iD[|V'^)(V'.|] = (3.38) 
spans a DFS if it is invariant under -ffcff- 



IV. GENERAL CONDITIONS FOR 

DYNAMICALLY STABLE 
DECOHERENCE-FREE SUBSPACES 

In the previous section we have studied pure, restricted 
DFS that fulfill the sufficiency criterion of Definition [U 
namely LY)[p[t)] = for all t. However, for unitary dy- 
namics, it is only necessary that the purity Tr[/9^(t)] of a 
state is preserved during the evolution (given an initially 



We show in Appendix [X] that all states statisfying Def- 
inition [2] undergo unitary evolution and therefore consti- 
tute a true DFS. 

Remark 5 . In Definition [2l it is absolutely crucial that 
the purity of the state p{t) is preserved for all times t, 
since non-unitary evolution is possible at intermediate 
times even when the initial state and final states are pure. 
For example, consider a 2-level atom subject to decay in 
the upper level by spontaneous emission. For this sys- 
tem, an initially pure excited state will eventually evolve 
into a pure ground state. Thus, initial and final states are 
both pure. However, for any intermediate time between 
t = and t —i' oo, the state of the system will be an inco- 
herent mixture of ground and excited states. Therefore, 
despite starting out as a pure excited state, this system 
will not remain pure at all times and will therefore not 
be classified as DFS by Definition [2l 

Using Definition[2]and Eq. (|2.4|) . it immediately follows 
that [23, Hi 



dtTv[p^it)]^2{Lu[p{t)]). 



(4.2) 



Therefore, DFS exist when (XdM) = 0, which is less re- 
strictive than LdIp] = 0. We will analyze the difference 
between these two conditions in Sec. [V] In the remain- 
der of this section we present a theorem that gives nec- 
essary and sufficient conditions for the existence of DFS 
according to Definition [2] and also provide a constructive 
protocol for determining both when DFS exist and the 
exact specification of all states within the DFS. 

Theorem 6. Let the time evolution of an open quantum 
system described in a finite- dimensional Hilbert space be 
governed by Eq. {2.4-^ with time-independent H^ff. The 
space 

span{\ipi), |?/>2), • ■ • , l^if)} = Ejurap (4.3) 
is a DFS satisfying Definitions^ iff the basis vectors fulfill 
Jil^k) = ci\i^k) (4.4) 
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for all I = 1, . . . , M and for all k = 1, . . . , K and Ej^mp 
is invariant under 



M 



(4.5) 



1=1 



This guarantees identification of all DPS as well as ex- 
plicit construction of all DF states. 

Proof. Before proceeding with the proof, wc point out 
that -Ejump contains only those eigenvectors of J; that 
span an invariant subspace with respect to H^v, i-e., the 
dimension of i?jump may be equal to or less than the 
number of simultaneous linearly independent eigenvec- 
tors of the jump operators Ji {1 = 1,..., M) and there is 
no requirement in the theorem that these be the same. 

We begin the proof by observing that for any M- 
dimensional vector b with complex coefficients hi, - ■ ■ , bu 

M 



. M 

1=1 



(4.6) 



where 



1=1 



J/(6z),pj/(6/)l + \jiibi)p,jl{bi) 



(4.7) 

with Ji{bi) ^ Ji- hi. 

Now we prove that the conditions of Theorem [5] are 
sufficient. Any state 



\i'{t)) e Ej^^p 



(4.8) 



can be written as 



\m) = ^aA:(i)l^fe> and pit) = m)){m\- (4.9) 



fc=i 



Then for c~ (ci, • • • , cm) 



Thus Ej^^p is a DFS. 

To prove that the conditions are necessary we suppose 
£^jump is a DFS. For fixed to, 

K 

l^(io)) =^afe(to)IV'fe), p(io) = |V'(to))(V'(<o)|, (4.13) 

k=l 



which implies 



9tTr [p2(io)] 

2Tv[p{to)LD[p{to)]] 
mo)\LD[p{to)mto)) 



(4.14) 



Suppressing the dependence on t^ we see that Eq. (|4.14[) 
has the same form as Eq. p.l4p . Repeating the same 
argument we are led to the condition = c/|V')- Fur- 
thermore, the line of reasoning presented in Eqs. p.28p 
and (|3.29[) proves that the eigenvalues c; are equal for all 
basis states IV'fc)- Thus we can conclude that j/|V'(to)) = 
Q(io)|V'(*o)), for any to- 

At this point, it appears that the eigenvalue ci{t) may 
depend on time t, but we argue now that ci{t) must 
be time-independent for finite-dimensional systems, i.e. 
ci{t) = Q. Eq. ()2.4p is a system of linear, first-order dif- 
ferential equations with constant coefficients; the com- 
ponents of p{t) and iV'(i)) are therefore continuous func- 
tions in time. Ji\ilj{t)) = c/(i)|'0(t)) then implies that 
ci{t) must be continuous, too. On the other hand, the 
time independent operator Ji acts on a finite-dimensional 
Hilbert space so that its spectrum is discrete. Thus, c;(t) 
must be constant. 

According to Eq. (|4T0)) for p{t) = \iit)){tp{t)\ with 
J;|i/'(i)) = ci\i]:{t)), we then have 

p{t) = -i[i?eff,p(i)] +iD[p(t)] = -i[H,,,p{t)] (4.15) 
and 



p{t) = exp(— i/fovi)p(0) exp(iifcvi) 



(4.16) 



M 



1=1 



. M 



1=1 



with iV-'(i)) = exp(— ii7ovi)|V'(0)) because Hcv is time- 
LdW)] = L^§'[p] - i[^ J2 ^' (^i*-^' " ^'^/) ' independent. Thus \iP{t)) e Ej^^p if Ej^^^p is invariant 

under Hcv The Hilbert space Tis can in principle con- 
tain multiple DFS. Each DFS is characterized by an M- 
tuple (ci , . . . , Cm) of eigenvalues of all the jump operators 
Ji, . . . , Jm for states in i?jump- 

Now we demonstrate that Theorem |6] is constructive, 
i.e., we explain how to use Theorem [6] to determine all 
basis states in all DFS. Condition 1 of Theorem [51 i.e., 
Eq. (|4.4p . states that DFS are subsets of all common de- 
generate eigenstates of the jump operators Ji, . . . , Jm- 
We label the space of all common degenerate eigenstates 
of the jump operators as Ej, with dimension N. The 
space Ej is determined by a collection of corresponding 
eigenvalues of the jump operators, ci, . . . ,cm- A differ- 
ent collection of eigenvalues might give rise to a different 
nonempty space Ej and constitute potentially a different 



because Ji{ci)\ilj{t)) = (J; - cii)\ip{t)) = 0|V'(t))- Hence 
the time evolution for p(t) is given by 



m 



Hcv,p{t) 



(4-11) 



so that 

9,Tr[p2(t)] = 2Tr[p(t)p(i)] 

= 2TrU[ifev,p(0]p(0 =0- (4.12) 



8 



DFS. In order to satisfy condition 2 of Theorem [6l i.e., 
invariance under Eq. (j4.5p . we need to restrict the space 
Ej to a subspace Ejump that is invariant under ffev. We 
accompHsh this task using Proposition O Proposition 
[5] says that a space that is spanned by some or aU of 
the collection of degenerate simultaneous eigenstates of 
a collection of operators ( Ji, . . . , Jm in this case), is in- 
variant under the action of a Hamiltonian H iff all states 
in the space are eigenstates of the corresponding com- 
mutator operators Bi = Ji] with eigenvalues 0, for all 
1 = 1,..., M . Thus, in order to construct DFS from Ej, 
we i) consider all basis elements \ipi), . . . , \'4'n) for Ej, ii) 
form linear combinations of these as ai • ■-\-aN\'4^N) 

with ai, . . . , Q!Ar free parameters, and iii) then determine 
all possible sets {a^, z = 1, . . . A^} such that the resulting 
states are eigenstates of Bi = [iJ, Ji] with eigenvalue 
for alH = 1, . . . , M . The resulting states form a complete 
basis for the DFS specified by eigenvalues ci , . . . , cm of 
the jump operators. Different choices of eigenvalues for 
the jump operators can give rise to different DFS and are 
constructed similarly. □ 

Remark 6. Theorem |6] applies only to finite-dimensional 
spaces. To demonstrate this we consider a damped har- 
monic oscillator described by the master equation 

p{t) = —i[u}{)a) a, p{t)\ + ^[2ap(i)a^ — a) ap{t) — p{t)a) d], 

(4.17) 

with a the annihilation operator. A coherent state 
pit) = \a{t)){a{t)l with a{t) = c-^"^"-''/^)* ^nd a 
a complex constant, is a solution to this master equa- 
tion. \oi{t)) is an eigenstate of a with unit purity, 
but it is not invariant under the evolution Hamiltonian 
7?ev(i) = t^oa^a + ma*{t)a - a{t)a)). Theorem i] does 
not apply to this example because a coherent state \a{t)) 
is an eigenstate of a with a time-dependent eigenvalue 
a{t). Thus, the analysis following Eq. (|4.14p in the proof 
of Theorem [5] does not apply. 

We now provide a constructive protocol for DFS, i.e., 
we present a procedure that explicitly determines the ex- 
istence and composition of all DFS. 

Protocol 1. 1. We consider all Lindhlad operators, 
Ji 1 = 1,..., M , from Eq. \3.5]) and calculate com- 
mon eigenstates for these operators, i.e., eigen- 
states satisfying Ji\ipi) = ci\ipi), fori = 
1,...,M. Note for each I, all \ipi) are eigenstates 
of Ji with the same eigenvalue. However, eigenval- 
ues for different index I need not be the same. The 
number of such common eigenstates is N > 0. For 
N = 0, there is no DFS. 

2. When N > we construct a space Ej = 
span{\il'i), i ~ 1, . . . , N}. This space is character- 
ized by a set of eigenvalues ci, . . . , cm of the jump 
operators. According to Theorem\^ any DFS must 
be a subspace of Ej . Note that according to Theo- 
rem\^ Dim{HDFs) < 



3. Now we apply Proposition [5| to the conditions de- 
termined by Theorem to select out the states in 
Ej that belong to a DFS (note that there is no re- 
striction to a single DFS). The conditions of The- 
orem require that the DFS be invariant under 
According to Proposition [31 a state \ipi) in 
Ej that is invariant under some Hamiltonian H is 
a zero eigenstate of the operator [H,Ji]. Thus, the 
DFS is a collection of simultaneous eigenstates of 
Ji that are also eigenstates of all Bi = [Hey,Ji], 
with eigenvalue 0, i.e., Hdfs = {IV'i) • G 
Ej and Bilipi) = Olipi)}. We can therefore sim- 
ply construct the operators Bi = [Hgy, Ji] and use 
these to select out the states \'>pi) in Ejump that also 
fulfill Bi\tpi) = 0\ipi) . Doing this for all I gives the 
basis for the DFS that is specified by eigenvalues 
Ci, . . . , Cm; of total dimension Dim{E jump) < N . 

This protocol must be repeated for all possible com- 
binations of common eigenvalues of Ji, J2, . . . , Jm- This 
will generate all existing DFS. These results provide sim- 
ple constructive criteria for the existence of DFS and for 
exact specification of all DFS states for any given Marko- 
vian master equation [4l| . 

V. INCOHERENTLY GENERATED 
COHERENCES 

In this section we analyze the difference between the 
two definitions of decoherence free subspaces introduced 
in this paper, i.e.. Definitions [1] and [H The following 
example demonstrates that there are states that satisfy 
Definition [2] but do not satisfy the criteria of Definition 

m 

Example 3. Consider a two-level system with basis 
states |1) and |0) in presence of a single decoherence op- 
erator J = (7_|_ + (7z, A = 2, and Hamiltonian H = CFy. 
Then the state = |1) has non-zero action of Lo on 
it, with Lo[\ip){'ip\] = 'i'i^y Ai') iM] = ~'^x 7^ , and 
yet it corresponds to a stationary solution of the master 
equation that fulfills {Lo[p\) = 0. 

In absence of the Hamiltonian H, the super-operator 
[ I V') ("0 1 ] would therefore evolve the system into a state 
that differs from \^p) and that would in general therefore 
be subsequently affected by decoherence. However in this 
example, H exactly cancels the effect of the decoherence 
super-operator at t = and hence for all times, ensuring 
that lip) remains stationary and decoherence- free. 

The reason why \ip) is stable in this example is that 
Ld[p] initially acts coherently on j?/;) and the coherences 
generated by this are canceled by the action of the uni- 
tary term —i[H, ]. It is easy to see that the initial 
coherence is generated between the state = |1) and 
the orthogonal state |0) . However this cancellation is spe- 
cific to these two states and does not hold for an arbitrary 
state of the two-level system, for which additional terms 
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of the form of Eq. (|4.7p will be present in the action of 
Ld[p]- Consequently, if the Hamiltonian term were ab- 
sent, Lu[p] would generate an incoherent evolution of the 
quantum state \^) because it only acts coherently on this 
particular state and its orthogonal complement, so that 
once the initial state has been infinitesimally coherently 
transformed to some superposition of these states, Ld [p] 
will gain some contribution of incoherent terms of the 
form of Eq. (|47)) . 

The existence of states that satisfy Definition [5] but 
not Definition [1] is surprising because it disagrees with 
the following, physically intuitive argument. Inserting 
the decoherence operator T into Eq. (j3.5p we can rewrite 
the master equation Eq. (|2.4p as 

M 

p = - iHnhP + ipH\ + ^\iJipj} (5.1) 

1=1 

with the non-Hermitian Hamiltonian H^^\^ = H^s — ir/2. 
In the context of quantum trajectory methods [29|, 
ifnh generates a continuous evolution in time, while 
^f^i ^iJipJi is interpreted as generating sudden quan- 
tum jumps at random times. The dynamics in absence 
of quantum jumps is given by 

p(t) ^ e-"^"V(0)e'*^n» . (5.2) 

Since F and HcS are both Hermitian operators, they cor- 
respond to the anti-Hermitian and Hermitian parts of 
Hnhi respectively. Hence, HcS cannot eliminate the effect 
of the decoherence operator F. Naively, one also would 
expect that the quantum jump term X]f=i ^iJipjj can- 
not cancel the effect of F. This argument leads to the 
conclusion that setting L-£,[p{t)] = V t is the only way 
to achieve unitary evolution. Example [3] demonstrates 
that this argument is not correct in general. 

The existence of these states can be quite generally 
explained with reference to the properties of Markovian 
master equations. Eq. (|4.6p shows that the Markovian 
master equation. E g. (12.41) . is invariant under the follow- 
ing transformation [3lj 

J, ^i, -6;/ for / = !,..., M (5.3) 

. M 

H,s^ H,s+^Y,Xi(btJi-biJiy (5.4) 
1=1 

Here bi ranges over all complex numbers. 

This transformation implies that Lolp] can always 
rewritten as a sum of two terms (see Eq. (|4.6p ). namely, 
i) a decoherence super-operator with respect to the 
transformed jump operators Ji — bil, and ii) a com- 
mutator between p and the Hermitian operator Ho = 
I i^i^i ~ ^ijf^ ■ Thus, the general decoherence 

super-operator [p] always contains a commutator with 
a Hermitian operator and, therefore, can in principle 
generate unitary evolution under suitable circumstances. 



Consequently, we cannot simply interpret i[_ffcff, p] and 
Ly)[p\ in the Markovian master equation, Eq. (|2.4p . as 
terms that always generate unitary and decohering dy- 
namics, respectively. 

Under special circumstances, in particular, when the 
state p is generated by a common eigenstate of all jump 
operators Ji, Eq. (|4.7p vanishes and the decoherence 
super-operator consists then of just the commutator with 
Hermitian Hq: 

LD[p\^-i[HD,p]- (5.5) 

In this situation, a state p decoheres because the Hermi- 
tian Hamiltonian Ho generated by Lo\p\ drives state p 
out of the subspace where its purity is preserved, and not 
because L u [p] causes a genuine decay in a quantum sys- 
tem. It is then possible that the effective system Hamil- 
tonian J?eff might compensate for such unitary leakage, 
resulting in preservation of the subspace. This is the situ- 
ation described in Example [3] above. For a given system, 
it may therefore be possible to construct or to modify 
the effective Hamiltonian specifically to ensure negation 
of the effects of Hq and to thereby stabilize states p 
in the subspace corresponding to one or more common 
cigenstates of the jump operators J; . 

Example [3] thus illustrates an intriguing general phe- 
nomenon, i.e., that even though decoherence in an open 
system S is caused by the interaction between system S 
and its reservoir, decoherence can be completely elimi- 
nated by a particular system Hamiltonian under special 
circumstances. We would like to emphasize here that 
this effect is very different from well-studied decoherence 
suppression phenomena such as dynamical decoupling 
[l7, 18„ 19„ i2ft] . In dynamical decoupling, decoherence 
effects are suppressed because the time scale for the sys- 
tem Hamiltonian is much faster than the timescale for the 
decay processes. The system Hamiltonian changes the 
state of the system so rapidly that decoherence processes 
are adiabatically eliminated. In contrast, in Example [3] 
the decoherence effects simply do not take place, as in 
any DFS. 

To understand more generally under which conditions 
the incoherent part Ld [p] of the master equation can sup- 
port a DFS through the generation of coherences, we in- 
troduce now the following special class of DFS states. 

Definition 3. States with incoherent generation of co- 
herences (IGC) are pure states p(t) that evolve unitarily, 
dtTr[p^{t)] = OVt > 0, but that have LD[p{t)] ^ at 
some time(s) t. 

Theorem [S] allows us to investigate the difference be- 
tween restricted DFS states in the sense of Definition [T] 
and subspaces composed of IGC states, which we shall re- 
fer to as IGC subspaces. Assume state p — fulfills 
Ji\iP) = qIV^) so that {Lu[p]) = dtTiipHt)] = 0. Then 

= E (2|cd'p - ciJIp - clpJi) . (5.6) 
The following result can be inferred from this equation. 
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Corollary 7. An IGC state can only exist if, for at 
least one I, Ji\tl') ~ qIV") with ci 7^ 0, and \^p) is not an 
eigenstate of J2i '^i^i ■ 

Remark 7. This result implies that IGC states and IGC 
subspaces are rare in nature. In particular, IGC states 
do not exist when the operators Ji are normal ([J;, jj] = 
dl]) since in this case any G Hufs would nec- 
essarily also be an eigenstate of j/. This implies that 
IGC states do not exist in systems subject to decohcrence 
due to dephasing (see also the discussion below Corollary 
3). In decoherence processes that involve population ex- 
change, such as spontaneous emission or thermal excita- 
tion, a DFS state typically corresponds to an eigenstate 
of Ji with Q = 0, as discussed earlier, and so Corollary[7] 
precludes any IGC in this situation. Thus IGC states do 
not exist in systems with decoherence due to indepen- 
dent decay, excitation or dephasing processes. Hence, 
while the intuitive physical idea that ioip] = should 
be fulfilled for DFS states is mathematically wrong, it 
nevertheless is valid for many physical applications. 

A. Examples of IGC states 

While IGC states do not exist for the most common 
physical decoherence models, they nevertheless can occur 
under special circumstances as the following examples 
illustrate. 

Example 4. Driven two-level atom in squeezed vacuum 
reservoir 

We consider a two-level atom that interacts with a ra- 
diation field. We also assume that the radiation field is 
prepared in a squeezed vacuum state. Then this system 
is described by the master equation (sij 

p = -t[Hs,p] (5.7) 
7o 2 / I o \ 

+ jo{sca-^pa^ + sc(T+pa+) , 

with s = sinh(r), c = cosh(r) and r the (real) squeezing 
parameter. Eq. (|5.7[) assumes Lindblad form when we 
introduce the operator 

J = ca^+sa+. (5.8) 

For r ^ the operator J has two non-zero eigenvalues 
±y^ with eigenstates \^p±) — ±-y/s|0) + -y/c|l). Each 
corresponds to a stationary IGC state if the Hamiltonian 
is given by Hs = ±^Vsc(s - c)ay, i.e., the DFS (IGC 
subspaces) are at most one-dimensional. 

However, one state is not sufficient to encode a qubit. 
Next we will explore how one can generate an IGC sub- 
space large enough to be used for quantum information 
processing. 



Example 5. N two-level atoms in the Dicke limit in 
squeezed vacuum reservoir 

We consider N two-level atoms in a squeezed vacuum 
reservoir in the Dicke limit [s^ l . We model the dynamics 
of this system by the master equation 

p = -i[Hd,p] + ^[^JpJ^ -J^Jp- pJ^J] (5.9) 

with 

N 

J ^'^Jn and J„ c(T„_ + so-„+, (5.10) 

11=1 

where the index n refers to the nth atom. 

To find IGC states, we need to construct all eigenstates 
of the operator J. They are given by all possible tensor 
products of iV'i)- For a given eigenstate, let n+ be num- 
ber of \'4>+) components and n_ = N — 71+ the number 
of IV^-) components. Then the corresponding eigenvalue 
is (n+ — n^)^/sc. For n+ = n_ the eigenvalue is zero, 
indicating that the state is not of IGC type. All other 
eigenstates of J become IGC states if the driving Hamil- 
tonian assumes the form = ■^(n+ — n_)y'ic(s — c)Sy 

with 4 = 

States with the same have the same eigenvalue and 
form an IGC subspace. The dimension of the IGC sub- 

f N\ 

space is given by the binomial coefficient: ( j. For 

> 3, there exist IGC subspaces with dimension greater 
than 1 and can be used to encode quantum information. 
For example, for n = 3 the IGC subspace correspond- 
ing to the eigenvalue ^/sc can be spanned from the non- 
orthogonal basis states: 

IV^i) = - V^(lioi) + |oii)) + \/^(|ioo) + |oio)) 

IV'2) = - cV^liii) + sV^lioo)) - cV^loii) + sV^lOOO)) 

K-g) =V^(|011) - |110)) - x/i(|001) - |100)). (5.11) 

The encoding efficiency of the IGC subspace (the number 
of logical qubits, log2(Dim(IGC(iV))y) approaches unity 
just as in the case of the usual DFS [3, 9- 

These examples seem to suggest that in finite- 
dimensional systems IGC states do only appear in rel- 
atively exotic situations. However, it is quite simple to 
find realistic examples in infinite systems. We give two 
examples below, noting that since Theorem 6 does not 
apply to infinite dimensional systems, these states may 
not be invariant under time evolution with iJ, despite 
the preservation of the state purity. 

Example 6. Driven damped harmonic oscillator 

Consider a driven harmonic oscillator with annihilation 
operator d and 

Ld[p\ = [2dpa^ — a^ap — pa^a] (5.12) 
iJ^g = ujQ{t) a) a + g* {t) a + g{t) a) . (5.13) 
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This model is frequently used to describe photons in a 
lossy single-mode cavity that is coherently pumped by a 
driving field of amplitude g{t). A coherent state |a) is an 
eigenstate of a with eigenvalue a and fulfills L/j [ja) (aj] =^ 
for a 7^ 0. In absence of a driving field {g{t) = 0) and 
for constant frequency ujo, this state remains a coherent 
state [3l|. This also holds in the general case; for a given 
function g(t) and initial condition a{0), a coherent state 
evolves as the pure state p{t) = \a{t)){a{t)\ with 

ait) = a(0)e^(*) - z f dt' e^^''>-^^'"> g{t') , (5.14) 
Jo 

and F{t) = J* dt' {—iujo{t') — 7/2). Hence, purity is pre- 
served and a coherent state corresponds to an IGC state. 

As we noted before, Theorem [6] does not apply to 
infinite-dimensional systems. Thus, purity is preserved 
for coherent states, although they are not invariant un- 
der the evolution Hamiltonian. 

For a coherent state, the effect of the decoherence term 
is reduced to an attenuation of the oscillator: it describes 
loss of energy (or photons) but not loss of coherence. 
However, if the choice oi g{t) is consistent with the trans- 
formation introduced in Eq. (j4.10p . for instance in the 
case u!o{t) = 1^0 and g{t) = ^jae'"^^"* , then purity and 
mean photon number (a^ap(t)) are both preserved, i.e., 
such a state does not experience decoherence or attenu- 
ation. 

Example 7. Two-photon absorber 

A two-photon absorber pumped by a two-photon 
parametric process constitutes another example of IGC 
states. The master equation for this system is 

p = -i[HD,p\ + ^(2aVat^ - o)'^a^p - po)'^a^) (5.15) 

with Hd = ^((a*)^a^ — a^a^^). The two coherent states 
Izba) arc solutions to this master equation [ssl. [s^. [bsI. [36| 
and arc both cigenstatcs of the two-photon annihilation 
operator a} with the same eigenvalue c? . Hence they 
form a two-dimensional IGC subspace. 



also provided a simple constructive protocol for explicit 
computation of states in all such DFS that exist. 

DFS states that fulfill (LD[p(t)]) = but for which 
^iAp{ty\ ^ are especially interesting. These states sat- 
isfy Definition [5] but not [T] We have shown that these 
states rely on incoherent generation of coherences (IGC). 
This means that they correspond to a set of states on 
which the decoherence term L£)\p(t)\ acts like a Hamil- 
tonian term, generating coherences with other states that 
are cancelled out by the system Hamiltonian. We showed 
that existence of these IGC states constitutes a general 
phenomenon that can be understood in terms of basic 
properties of the Markovian master equation. In par- 
ticular, their existence was shown to result from an in- 
variancc of the master equation with respect to a con- 
tinuous transformation, implying a dynamical symmetry 
that links the coherent and decoherent terms and that 
results in cancellation of the decoherent terms in certain 
circumstances. 

IGC states form a subset of DFS states and as such 
may be useful for quantum information processing. A 
general strategy to generate such states is to first identify 
states on which the decoherence superoperator acts co- 
herently for infinitesimal times and then to stabilize these 
states to make IGC states by adding a Hamiltonian to the 
master equation that cancels the action of Ljj [/^(t)] on the 
states. We have studied under which decoherence condi- 
tions a system admits IGC states. For finite-dimensional 
systems the most common decoherence models, includ- 
ing independent dephasing and incoherent excitation or 
de-excitation processes, do not allow for IGC states. In 
other words, the set of restricted DFS is then identical 
with the full DFS. IGC states in finite-dimensional sys- 
tems can only appear in relatively exotic situations, for 
instance when the system interacts with a reservoir that 
is prepared in a squeezed vaccum state. On the other 
hand, in infinite-dimensional systems it is easy to find 
examples of IGC states, e.g., coherent states in a lossy 
single-mode optical cavity. 



VI. SUMMARY AND DISCUSSION 

We have examined two different definitions of states 
p(t) of a quantum system that undergo unitary evolution 
in the presence of decoherence effects, both of which are 
related to the preservation of state purity. For the re- 
stricted DFS introduced in definition [1] the decoherence 
term in the master equation (j2.4[) vanishes, LY)\p{1b)\ = 0. 
This is a sufficient but not a necessary criterion for uni- 
tary evolution. General DFS states can be characterized 
by the less stringent condition, dt^T\p^[t)\ = for all t > 
0, or equivalently, {LY)[p{t)]) = 0, that was introduced in 
definition [21 With theorems H] and [H] we have established 
rigorous conditions for the existence of dynamically sta- 
ble DFS of both types in finite-dimensional systems. We 
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APPENDIX A: EQUIVALENCE BETWEEN 
UNITARY DYNAMICS AND PURITY 
PRESERVATION 

Here we show that Definition [2j according to which 
a state is DF if its purity does not change from initial 
unit value, is equivalent to requiring unitary dynamics for 
Markovian master equations. We start with the fact that 
the evolution of any quantum state po is always described 
by a completely positve trace-preserving linear map (CP 
map) Stipo) ■ For such a map the following proposition 
holds. 

Proposition 8. Let pQ be a density matrix with purity 
one, i.e., Tri^p^) = 1. Then the action of a linear CP 
map Stipo) is unitary, 

£tiPo) - UtPoUl , (Al) 
ijfStipo) preserves purity, Tr(£t{po)'^ )=1. 

Proof. We first show that the condition is sufficient. As- 
sume that Stipo) = UtPQU}. Then 

TY{£t{pof) = TviUtPoUjUtPoUj) (A2) 

= Tr(C/tpoPof//) (A3) 

= Tr(poPo) (A4) 

= 1 . (A5) 

To prove that the condition is necessary, suppose that 
Tr(£'t(yOo)^) = 1 for a given pure state po- £t(pp) can be 
written in the operator sum representation as [38| 

£t{po) = J2^^po^l ' E^I^^- = i- (A6) 

k k 

For a state of purity one, the CP map on the pure state 
Po = can then be written as 

£mm = Y.m{M , (A7) 

k 

This is a uniform mixture of states \i}}k)- Such a mixture 
can only be of purity one if all \il)k) are proportional to 
each other, i.e., if 

l^fe) = A,#), (A8) 

where is a normalized state defined as = 
Ei\i}})/\\Ei\i}j)\\, provided that \\Ei\ip)\\ ^ 0. If this is 



not the case, the normalized state can be defined in terms 
of any other Ek for which 1-0)11 ^ 0. There must be 
at least one such operator E^, since if there is no such 
operator, then £t{po) = 0, which would contradict the 
assumption that purity is preserved. Using Eq. (jASp in 
Eq. (|A7|) . we then obtain 

^mW) = 5]£;;c|0)(V'|£;l (A9) 

k 

= Y.\^k\'\i')m (AlO) 
Now we use the unit normalized state \ip) to write 

Ei^^i' = E(^'i^fc^^i^) (All) 

k k 

= Y^i^PlElEkW (A12) 

k 

= im^') (A13) 

= 1, (AM) 
so that Eq. (jAlOp becomes 

Stmiijl) = |0)(0|. (A15) 

Hence £t maps \ip) to {ip). Such a map between unit 
normalized pure states can always be written as a unitary 
transformation Ut and explicit construction of Ut can be 
made, e.g., by using a Householder reflection [3^. Hence 

f*(|0)(0|) = [/,|^)(0|C// . (A16) 

□ 

Under the usual assumption of and initially factor- 
izable system and environment state, the time evolu- 
tion of a system state p{t) can always be described by 
a time dependent CP map £t{po) acting on the ini- 
tial state po [111 . By proposition [51 requiring unitary 
dynamics for all t is then equivalent to the condition 
Tr(p2(i)) = 1V< > 0. For dynamics that are governed 
by a Markovian master equation (|2.4p . p{t) is a diff'eren- 
tiable operator-valued function. Hence the time deriva- 
tive of p(t) is well-defined and purity preservation can 
be expressed as dtTi{p'^{t)) = OVt > with Tr(/9§) = 1. 
Therefore, for Markovian dynamics. Definition [2] is equiv- 
alent to requiring a DF state to undergo unitary dynam- 
ics. 
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